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. Abstract 

We explore in some detail the scenario proposed to explain the observed 
knee of the cosmic ray (CR) spectrum as due to the effects of photo disintegra- 
. tion of the CR nuclei by interactions with optical and soft UV photons in the 

' source region. We show that the photon column densities needed to explain 



' the experimental data are significantly lower than those obtained in previous 

estimations which neglected multinucleon emission in the photodisintegration 
O I process. We also treat more accurately the photodisintegration thresholds, we 

^ I discuss the effects of photopion production processes and the neutron escape 

■ mechanism, identifying the physical processes responsible for the qualitative 

features of the results. This scenario would require the CR nuclei to traverse 
\ column densities of ~ 5 x 10^^ — 2 x 10^^ eV/cm^ after being accelerated in 

' order to reproduce the observed knee, and predicts that the CR composition 

^ should become lighter above ~ 10^^ eV. 



1 Introduction 

It is widely accepted nowadays that cosmic rays (CRs) with energies per particle 
up to about 10^^ eV are protons and nuclei of galactic origin. Furthermore, it is 
well established that the full CR energy spectrum has a power-law behavior with a 
steepening taking place at the so-called knee, corresponding to an energy Eknee = 
3 X 10^^ eV. Although it is well known that the CR composition below the knee 
has a significant heavy component, its behavior beyond the knee remains somewhat 
controversial. Indeed, above 10^^ eV/nucleus the information about the CR mass 
composition has to be drawn from extensive air shower observations at ground level 
with the results from different experiments not always being compatible. Moreover, 
the theoretical predictions depend on the hadronic model adopted and this can 
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introduce further uncertainties on the inferred composition. The model dependence 
in certain observables can however be useful to discriminate among the hadronic 
models using the experimental data, and in this comparison QGSJET and VENUS 
turn indeed to be favoured against other models such as SIBYLL Thus, the 
CR mass composition beyond the knee is not definitively established yet, with some 
observations [0] suggesting that it turns lighter and others instead suggesting that 
the heavier components become dominant |, |], H- 

The attempts made so far to explain the physical origin of the knee of the cosmic 
ray spectrum can be roughly classified into three kind of models. One of them ex- 
ploits the possibility that the acceleration mechanisms could be less effective above 
the knee |0, ^, , while a second one assumes that leakage from the Galaxy plays the 
dominant role in suppressing CRs above the knee |]TU|, |TT|, |T2|. The third scenario, 
originally proposed by Hillas [|13|, considers nuclear photodisintegration processes 
(and proton energy losses by photomeson production), in the presence of a back- 
ground of optical and soft UV photons in the source region, as the main responsible 
for the change in the CR spectrum. The first two scenarios are based on a rigidity 
dependent effect which consequently produces a change in the spectral slope of each 
nuclear component with charge Z at an energy ~ EkneeZ. Hence, it removes first 
the protons and the lighter nuclei and only at larger energies affects the heavier 
nuclei, predicting that the CR mass composition should become heavier beyond the 
knee. On the contrary, the third scenario predicts a lighter CR composition above 
the knee due to the disintegration of the propagating nuclei. The predictions of this 
last scenario have been worked out by Karakula and Tkaczyk some years ago. 
The aim of the present paper is to re-evaluate in detail the propagation of nuclei 
and protons in the source region following their approach, but taking into account 
some additional features that turn out to be important for the determination of the 
source parameters required for the scenario to work. In particular, we re-evaluate 
the nuclear photodisintegration rates taking into account multinucleon emission pro- 
cesses, we introduce more accurate threshold energies in the giant dipole resonance 
for every stable nucleus (according to recent remarks by Stecker and Salamon [0), 
and discuss the impact of the neutron mechanism |]16[ which would allow the escape 
of neutrons from the source without further energy losses. Moreover, we show that 
by introducing a non negligible lower cutoff to the power-law photon distribution, 
the abrupt suppression of the all-particle CR spectrum above 10^'' eV previously 
obtained |]T^ is prevented, and hence the underlying reason for that effect is iden- 
tified. Let us also mention that although we focus our study here on the scenario 
proposed to explain the observed steepening at the knee, we expect that our general 
considerations about the treatment of photodisintegration processes should also be 
relevant in other contexts. 
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2 The propagation of cosmic rays 



2.1 The source and the photon background spectra 

As already mentioned, we deal with a model that considers that the CR nuclei are 
accelerated inside discrete sources which are surrounded by a very strong background 
of optical and soft UV photons. If the typical energy of the photons is in the optical 
range (1 — 10 eV), the photodisintegration of CR nuclei will start to be efficient 
at CR energies E > A x 10^^ eV. This is so because in the CR rest frame the 
optical photon (which is boosted by a relativistic factor 7 > 10^) appears as an 
energetic gamma ray with £^ > 1 — 10 MeV, i.e. capable of photodisintegrating the 
nucleus. Hence, the basic assumption is that the CR spectrum emitted by the source 
is a featureless extrapolation of the spectra measured below the knee and that the 
observed change in slope is actually due to the CR interactions (i.e. the process 
of nuclear photodisintegration and the photopion production by protons) with the 
surrounding photon background. 

We then assume that the source emits nuclei with mass number A (with A 
ranging from 1 to 56). Since there is only one stable isotope for a nucleus of a given 
mass A along most of the decay chain from ^^Fe to ^H, we take for definiteness a 
unique charge Z associated to any given value of A. The differential fluxes emitted 
by the source are assumed to be given by power-law distributions 

4>^{E) = (1) 

(where i = 56 — A throughout). The intensities $° and spectral indices 7^ were 
taken from the detailed knowledge about CR mass composition below the knee (for 
energies per particle above ~ few Z x 10^° eV) [0. 

For the photon background that surrounds the source we consider that its spec- 
trum follows either a (thermal) Planckian-type distribution given by 



' ' 2C(3)(A;BT)3exp(e/A;Br)-l 

where ks is the Boltzmann constant, T the absolute temperature and ( the Riemann 
zeta function (such that ({3) ~ 1.202), or a power-law distribution given by 

nie) = Pia - 1) (-^ - \- (3) 

where e^, are the lower and upper energy cutoffs and a the assumed spectral 
index. For both distributions, p represents the total number density and n{e) the 
differential number density corresponding to a photon energy e. Defining pe the 
total energy density, one has for the Planckian spectrum 
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PE^2mikBTp , (4) 
while the analogous expression for the power-law spectrum is 



We adopt a = 1.3 as in ref. |TJ], but the introduction of a non negligible lower 
cutoff em in the power-law spectrum excludes the presence of an abundant number 
of low energy (infrared) photons that would otherwise play the dominant role in the 
photodisintegration above 10^^ eV, as was the case in ref. [jl^. Regarding the spatial 



dependence of the radiation background we do not make any particular assumption, 
but the results will depend only on the integrated photon column density along the 
trajectory of the CR since its production at the source until it leaves the region filled 
with radiation. 



2.2 Propagation of nuclei: photodisintegration 

The main mechanism of energy loss for nuclei with energy per particle below 10^^ eV 
propagating through these photon backgrounds is the process of photodisintegration. 
In fact, a nucleus of mass A = 56 — i and Lorentz factor 7 = E/AnipC^ that 
propagates through a photon distribution n(e) (as, for example, those of eqs. (2) 
and (3)) has a probability of fission with emission of j nucleons given by 



^^.(^) = ^ /, d6 4^ de'6V,,(6'), (6) 



Je' ../27 

where aij is the corresponding photodisintegration cross section, e the photon energy 
in the observer's system and e' its energy in the rest frame of the nucleus. In order 
to calculate Rij we fitted a-ij with the parameters given in Tables I and II of ref. 
T8[] , while the reaction thresholds were taken from Table I of ref. [|1^]. A 



useful quantity to take into account all reaction channels is the effective emission 
rate given by 

Ri,eff = ^3Rij ■ (7) 

At low energies {e[j^^^ < e' < 30 MeV), the cross section (Tij{e') is dominated 
by the giant dipole resonance and the photodisintegration proceeds chiefly by the 
emission of one or two nucleons. Note that, according to recent remarks pointed 



out in ref. ||T5[, we employ threshold energies that depend both on the nucleus 
and on the number of emitted nucleons. The improvement consists in shifting the 
threshold energy, that was previously assumed to be e[f^j. = 2 MeV for all reaction 
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Figure 1: Plots of Ro^eff/PE versus energy per nucleon E* for a thermal photon 
spectrum (with ksT = 10 eV). -Ro^e// "^^^ evaluated using the fixed threshold energy 
= 2 MeV, while -Ro,e// corresponds to the improved threshold energies according 



to ref. \15 . 



channels and all nuclei |14, 1^, to higher energies, such that single- nucleon emission 
has now a typical threshold of ~ 10 MeV, while the double-nucleon emission energy 
threshold becomes typically ~ 20 MeV. Although the appearance of the knee in 
this scenario is ultimately a threshold effect, the improved calculation including the 
change in the specific threshold energies has however a minor impact on the resulting 
emission rates when compared with those obtained with a fixed threshold energy 
e[i^^ = 2 MeV. This can be appreciated in Figure 1, where plots of Ri^eff/pE vs energy 
per nucleon are shown for a ^^Fe nucleus {i = 0) propagating through a thermal 
photon spectrum (with ksT = 10 eV). The emission rate Ro,eff (calculated with the 
specific threshold energies) shows a small shift towards higher energies with respect 
to that obtained by means of the fixed threshold energy 

(labeled The reason 

for this reduced effect of the threshold shifts is that the dominant photodisintegration 
effects are due anyhow to the energies around the giant resonance peak, which for 
single-nucleon emission are typically peaked around ~ 20 MeV and for double- 
nucleon emission around ~ 26 MeV, and hence clearly well above the threshold 
energies. 

At higher energies, cTjj(e') is approximately flat and the multinucleon emission 
acquires a higher probability, becoming actually dominant for heavy nuclei, for which 
the probability for single nucleon emission is of only 10%. To estimate the relevance 
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of multinucleon emission processes, we compare the effective rate from eqs. (6) and 
(7) with the emission rate that neglects multinucleon processes (as calculated in 
ref. JT^)- Figure 2(a) shows plots of Ri^effl Pe versus energy per nucleon E* for a 
^^Fe nucleus {i = 0) propagating through a thermal photon spectrum (with ksT = 
1.8 eV). From the figure it becomes evident that multinucleon processes play an 
important role, increasing the emission rate at high energies by a factor of ~ 4 
with respect to the single nucleon emission results. Analogously, Figure 2(b) shows 
the corresponding plots for a power-law photon spectrum with upper energy cutoff 
eM = 20 eV and spectral index a = 1.3. From this figure we can also see the 
effect of introducing a non negligible lower cutoff in the spectrum. In fact, we 
observe that with a negligible cutoff (e.g. = 10~^ eV) the emission rate grows 
steadily, because increasingly abundant low energy (infrared) photons give in this 
case the dominant contribution to the rates. They also keep the influence of the 
giant resonance dominant even at high particle energies, and hence in this case the 
multinucleon emission has little impact. However, if we take a non negligible cutoff 
into account (for instance, = 0.8 eV), we find that the fiat high energy regime of 
(Tj now dominates beyond a given value of E* and hence the emission rate saturates 
at high energies. Comparing the two emission rate curves with non negligible cutoff, 
we observe again that multinucleon emission processes lead to a significant increase 
(by the same factor of ~ 4 we already encountered) in the emission rate. Thus, 
introducing a non negligible lower cutoff energy in the power-law photon spectrum 
the results should not differ much from those with the thermal photon spectrum 
and, being the rates with multinucleon emission larger, they should also require 



smaller photon densities than those previously estimated in ref. [14 to produce 
similar overall effects. 

It should be noticed that in Figures 1 and 2 we have chosen to plot the emission 
rates as functions of the energy per nucleon E*, instead of energy per nucleus E, 
because it is E* that remains constant during photodisintegration and also because 
for different nuclei the shapes of the cross section are similar with the maxima 
occurring approximately at the same values of E*. 

The equations that describe the propagation of nuclei are: 



= -UE*,x)J2R,,{E*) + {l-5oi)Y.Ri^-J)J{E*)ct>,^,{E\x), (8) 



where (f)i{E*, x) is the differential flux of a nucleus of mass A = 56—i (for < z < 54) 
with an energy per nucleon E* at propagation distance x from the source, Rij{E*) 
is the above defined emission rate and Sij is the Kronecker delta. The exact solution 
of these equations (i.e., the differential fluxes at propagation distance L) is given by: 
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Figure 2: Plots of Ro^eff/PE versus energy per nucleon E*. (a) Results for a thermal 
photon spectrum with ksT = 1.8 eV. (b) Results for a power-law photon spectrum 
with upper energy cutoff sm — 20 eV and spectral index a — 1.3. The multinucleon 
emission rate corresponds to a lower energy cutoff = 0.8 eV, while the single- 
nucleon emission rates correspond to = 10~^, 0.8 eV, as indicated. 
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where 



<l>i{E% L) = Y, bij{E*) exp(- ^ Rjk{E*)L), (9) 

j=0 k>l 



bu = (PKE*)-{l-So^)Y.b^AE*) (10) 

j=0 



and 



J2k>lRikiE*) - J2k>lRjkiE*) 

Recalling that the emission rates Rij are linear in the photon energy density 
Pe, it should be noted that, as expected, the fluxes only depend on the column 
density PeL. It also should be pointed out that, since the equations describing the 
CR propagation only depend on the energies per nucleon E*, the final CR spectra 
require the re-calculation of the solution (eqs.(9)-(ll)) in terms of energy per nucleus. 

2.3 Propagation of protons 

The main energy loss mechanism that must be taken into account for the propagation 
of protons is the photomeson production process. We have to distinguish between 
the primary protons produced at the source with an initial flux ^^=55 (as given by 
eq.(l)) and the nucleons released as by-products of the nuclear photodisintegration. 
The fluxes corresponding to these (secondary) protons and neutrons (0*^ and 0", 
respectively) can be easily determined from the solution for the nuclear fluxes, as we 
will see in the following. Let ct)f{E,x) be the differential flux of nucleons of energy 
E produced within a distance x away from the source by the photodisintegration of 
nuclei of mass number A — 56 — i and energy per nucleon E* — E/A. Then, the 
equation governing the evolution of the fluxes is 

^^l^l^^^^{E*,x)Y^jR,j{E*) , for0<i<54 , (12) 
ox .^^ 

(with the initial condition 4>f{E,x = 0) = 0) and its solution reads 



j=0 



' Y.k>lkR^k{E*) 

, Eik>i Rjk{E*) 



k,{E* 



l-exp[-J2R,k{E*)L 

k>l 



(13) 



Thus, the total fluxes of secondary protons and neutrons are found by performing 
the weighted sums over index i: 
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54 y 

0^nE,L)=^— ^0f(E,L) , (14) 

54 . 7 \ 

0'^(E,L)=5:(l-^)0f(^,^) , (15) 

where is the charge of the i-nucleus. 

We will also discuss the possible impact of the so-called neutron mechanism [|16||, 



which could take place in the presence of a magnetic field in the source region. 
Indeed, a magnetic field would retain the charged particles within that region for 
some period of time, while non-charged particles would be able to escape rapidly. 
In that case, we will consider that the flux of neutrons 0" adds to the final all- 
particle CR flux neglecting any energy losses, while the flux of secondary protons 
0*^, combined with the initial flux of primary protons 0^=55, will undergo energy 
losses due to the well studied process of photopion production. Indeed, since a 
background photon looks like a high-energy gamma ray in the proton rest frame, 
the photopion reactions 

p + 7 — s> J9 + 7r° (16) 
p + 7 — s> n + vr"'' (17) 
proceed with an interaction probability given by 

g{E) = — I de^ I de'e'a,,{e')K{e'), (18) 



where, analogously to eq. (6), 7 = E/mpC^ is the Lorentz factor of the proton, a^p 
is the 7p interaction cross section (parametrized from experimental data compiled 
in ref.|]19[), K is the coefficient of inelasticity (defined as the average relative energy 
loss of the proton), e the photon energy in the observer's system and e' its energy 
in the rest frame of the proton. The energy threshold for photopion production is 
e;,, = 145 MeV. 

Then, if we call (j)^{E,x) and 0*"(£',x) the differential fluxes of protons and 
neutrons of energy E (produced by the reactions given by eqs.(16) and (17), respec- 
tively) after tarversing a distance x from the source, the corresponding propagation 
equations are 

- .,^E)nE, .) + (r^) ^ A (19) 



dx V , . 2{1- K)'' \1- KJ ^ \1-K 

and (neglecting neutron energy losses) 

' ' Vf^.x) . (20) 



dx 2(1-/0 Vl-iry \l-K' 
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where the factors | arise from the relative probabihty of occurrence of reactions 
(16) and (17). The coefficient of inelasticity K was given a typical value K = 0.3, 
irrespective of proton energy |T^. The initial conditions to solve these differential 



equations are 

^P{E, x = 0) = <PU,{E) + 0^^(E, L) (21) 



and 



XE,x = 0)=0 , (22) 



respectively. In the first one we assumed for simplicity that the secondary protons 
produced by photodisintegration of nuclei add directly to the initial proton spec- 
trum. This is justified since whenever the photopion processes are relevant, the 
photodisintegration is very efficient and hence nuclei are rapidly disintegrated. The 
solution to eq.(19) turns out to be 

oo 

ctf{E,L) = Y.<f,{E,L) , (23) 

1=0 

where 

ct>l{E,L)=ct>nE,Q)eM-9{E)L) , (24) 

while for / > 

^iE L) - '^'(tA'-O) a „ / E \ ' exp (-9 (^) l) 
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(25) 



For the neutron component, the solution to eq.(20) reads 

00 

<P*-{E,x) = J2<Pr{E,x) , (26) 
1=0 

where 

while for / > 



L*n 



(2(i-ir))'+i f}.^ [{1 - Ky^' ' 
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I 

The linearity of the photopion reaction rate (?(-£') with respect to the photon 
energy density pe (as imphed by eq.(18)) causes the proton and neutron fluxes to 
depend only on the column density PeL, similarly as for the fluxes of nuclei. 



exp [-9 ( (i_^)n+i ) ^) 



(28) 



3 Results and comparison with observations 

So far, we wrote down the appropriate equations governing the propagation of CR 
particles, and we found their exact solutions: eqs.(9)-(ll) for nuclei, eqs.(13),(15) 
and (26)- (28) for neutrons (that eventually decay into protons outside the source 
region), and eqs.(23)-(25) for protons (that undergo energy losses due to photopion 
production), the sum of all these contributions being the total CR flux 4>totai- We 
can now investigate how 4>totai depends on the assumed photon distributions around 
the source, namely on the column density peL and on the parameters involved in 
the photon spectral distributions (eqs.(2) and (3)). 

Figure 3(a) shows the total CR differential flux as a function of the energy per 
particle for a Planckian spectrum with = 1.8 eV and several values for the 
column density peL. Analogously, Figure 3(b) shows the results corresponding to a 
Planckian spectrum with fc^T = 10 eV. To compare our results with experimental 
data, the figures also show the observed spectra measured by different experiments 
(CASABLANCA, DICE, Tibet, PROTON satellite and AKENO) 0, |, |^ . 



The effect of increasing the temperature is clearly that of shifting the knee to- 
wards lower energies, while the effect of increasing peL is that of intensifying the 
steepening at the knee. This behavior is easily explained by the fact that it is the 
nuclear photodisintegration the responsible for the occurrence of the knee (since 
photopion production processes show up at higher energies). Then, if the photon 
temperature (and thus, the mean photon energy) is larger, the particle energy re- 
quired for the nuclei to disintegrate decreases. Furthermore, CR propagation cannot 
be affected by peL below the knee, where the photodisintegration has not set up yet, 
but above the knee the disintegration rates have a linear dependence on pe- Hence, 
we should expect peL to control the steepening in the knee region, in agreement 
with the results shown in the figures. 

In ref. [|14| the best fit to observations for the Planckian spectrum was that 



corresponding to ksT = 1.8 eV and peL = 2.25 x 10^^ eV/cm^. For the sake of 
comparison. Figure 3(a) shows our results for the same temperature. It is clearly 
observed that the new results are compatible with the AKENO experimental data 
for low column density values peL ~ 5 x 10^'^ eV/ cm^, i.e. ~ 50 times lower than the 
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Figure 3: Plots showing the total CR flux versus the energy per particle for thermal 
spectra with (a) ksT = 1.8 eV and (b) ksT = 10 eV, and several values for the 
column density PeL, as indicated. For the sake of comparison, experimental data 
measured by different experiments are also shown. 
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previous result. Figure 3(b) shows that the corresponding results for fc^T = 10 eV 
are in somewhat better agreement with CASABLANCA observations for column 
densities in the range peL = 5 x 10^^ — 2 x 10^^ eV/cm^. One has to keep in mind 
however that the model discussed is oversimplified in many respects: we assume 
that the same photon column density is traversed by all nuclear species and for all 
CR energies, we neglect additional effects due to change in diffusion properties at 
these energies, which in alternative scenarios |12| are actually the responsible for the 
appearance of the knee, we assume that CR interactions in the interstellar medium 
are negligible, etc.. Hence, an accurate agreement with observational data is not to 
be expected, although it is reassuring to see that the results have the right tendency 
to account for the main qualitative features of the data in the knee region. 

Figures 4(a) and 4(b) show plots of mean mass composition (InA) versus E 
corresponding to the Planckian photon spectra with ksT = 1.8 eV and ksT = 10 eV, 
respectively. For comparison, we also show the CR mass composition measured by 
different experiments 0, ||, One can observe that solutions corresponding to 



the largest plotted values of column density peL imply that CRs above the knee 
are chiefly constituted by protons, in disagreement with experimental data; the 
non-negligible contribution of heavier CR components beyond the knee imposes an 
upper limit on the column density around {pEL)max ~ 2 x 10^^ eV/cm^. This result 
stresses again the fact that multinucleon emission processes play a significant role in 
CR propagation within this scenario. We also see that, although the solutions with 
lower column density {peL < {pEL)max) are in reasonable agreement with some 
data sets within error bars, the composition beyond the knee seems to be somewhat 
suppressed. Nevertheless, we should bear in mind that we have employed just one 
source with a fixed column density, and that in a realistic situation there would 
be many sources, each with a different environment of radiation. Avoiding any 
detailed calculations, we may think about the contribution of one additional source. 
If it happens to have a column density peL somewhat smaller, it would lead to a 
flatter flux; then, whereas the first source would be responsible for the occurrence 
of the knee, the second one would keep a non-vanishing nuclear component beyond 
the knee. Furthermore, in the computations we assumed that all nuclei traverse 
the same distance L, thus disregarding the fact that the magnetic field present in 
the source region should retain more effectively the heavier nuclei at any given CR 
energy. Taking this into account (for instance, assuming that a nucleus traverses a 
distance in the photon background proportional to its charge) we would find that, 
while the disintegration would be very effective for heavy nuclei (and would be, 
in fact, responsible for forming the knee), it would be less effective for the lighter 
nuclei (since the column density for them would be much smaller) and they could 
then survive beyond the knee. Hence, we see that this scenario can easily account for 
a change in the composition towards lighter nuclei above the knee, which is apparent 
in some of the existing data. On the contrary, rigidity dependent scenarios predict 
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the opposite trend, and may have difficulties to account for some measurements. It 
has indeed been pointed out that they would require the introduction of an ad-hoc 
additional light component above the knee to reproduce the observations . 

The results for the power-law photon spectrum (with a lower cutoff energy in 
the infrared range) should not differ much from the ones obtained for the Planckian 
spectrum, and this behavior was in fact corroborated. Let us fix the spectral index 



and the upper cutoff energy at the values found in ref. [|IJ] to correspond to the 
best ffi to experimental data, namely a = 1.3 and eA/ = 20 eV. Figure 5 shows (ptotai 
versus E for power-law distributions with column density peL = 10^* eV/cm^ and 
several values for the lower cutoff energy (em = 0.01,0.1,0.8 eV). It is clear that 
the overall spectrum displays a steepening which has the right qualitative features 
to explain the knee although it does not provide an accurate ffi to neither AKENO 
nor CASABLANCA data individually. Turning now to discuss the dependence of 
the results on the lower cutoff energy, it is found that, apart from minor differ- 
ences arising from spectrum normalization, lower values of yield lower fluxes in 
the high energy end. The explanation for this is that from the expression for the 
photodisintegration emission rate (eq.(6)) it turns out that the role of is that of 
determining a value of the particle energy E, such that below E the giant dipole 
resonance dominates, while above E the flat high energy regime of a prevails. As 
we lower e^, more photons with lower energy become available so that the flux 
suppression imposed by the giant resonance extends to higher particle energies. 

This analysis provides a suitable explanation of the results shown in Figure 5, as 
well as those of Figure 2(b). The results previously obtained for power-law spectra 



show an extremely abrupt flux suppression that seems to be in contrast with 



experimental data (CR energy spectra and, especially, CR mass composition) above 
E = 10^^ eV, and can be attributed to the negligible value for adopted there. 
The use of a non negligible cutoff prevents the flux suppression and provides a better 
fit to observations in this energy region. It should also be remarked that a solution 
yielding an excessive CR flux at energies above 10^^ eV is by no means troublesome 
since it is clear that at these high energies efficient leakage of CRs from the Galaxy 
cannot be disregarded, and that a less effective acceleration could also be playing 
a role. These mechanisms would clearly contribute to lower the CR flux, so that it 
does not seem reasonable to expect a reliable explanation of the CR spectrum by 
photodisintegration processes alone beyond E ~ 10^^ eV. 

Figure 6 shows plots of the proton and neutron components (0^-|-0^^ and 
respectively) as well as the total CR flux (ptotai for the Planckian photon spectrum 
with fc^T = 1.8 eV and peL = 2 x 10^^ eV/cm^. In order to estimate the possible 
effects of the neutron mechanism, we also plotted the total CR flux that is obtained 
ignoring the escape of neutrons, labelled (plotai- From the figure it is seen that the 
neutron mechanism raises the CR total flux by ~ 100% for ii^ > 5 x 10^^ eV. 
Had we adopted a column density PeL = 10^^ eV/cm^, the corresponding flux 
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Figure 4: Plots of mean mass composition (InA) versus E corresponding to the 
Planckian photon spectra with (a) ksT — 1.8 eV and (b) ksT — 10 eV. For com- 
parison, the CR mass composition measured by different experiments is also shown. 
Notice that several data sets exhibit the same qualitative behavior, since they cor- 
respond to the same experimental data and differ only in the hadronic models used 
to interpret the observations. 
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Figure 5: Plots showing the total CR flux versus the energy per particle for power- 
law spectra with a. — 1.3, eu — 20 eV, p^L = 10^^ eV/cm^ and several values for 
the lower cutoff energy (e^ = 0.01,0.1,0.8 eV). For comparison, experimental data 
measured by different experiments are also shown. 

increment would be only ~ 10%. Nevertheless, it should be noted that these results 
correspond to very high column densities (actually inconsistent with the observed 
CR composition) , and that effects arising from the neutron escape mechanism turn 
out to be negligible for lower values of the column density p^L, for which also the 
photopion losses of protons can be altogether ignored. 

To sec what arc the implications for the sources in order to produce the required 
photon densities for this scenario to work, let us make the following simple estimate. 
An isotropic source with photon energy density pe,o at its surface (with radius ro) 
will have around it a photon energy density given by pEif) = PE,o{ro/ry. This 
means that the photon column density measured radially from the surface is (pEf^) — 
PEfif^o- On the other hand, if the CRs are magnetically confined near the source, the 
path they travel across the dense photon field will be larger, i.e. (peL) = X{pEr), 
with the parameter A ^ 1. If to obtain an estimate of the maximum achievable 
densities we assume that the source radiates with a luminosity L = ATrr^pEfiC smaller 
than the Eddington limit, i.e. L < Lem = AnGMnipc/aT ~ 1.3 x 10^^M/Mq erg/s, 
with M the mass of the source, nip the proton mass and Gt Thomson's cross section, 
we find that the radius of the source needs to satisfy 

'^0 < -\t7{\\ 1 — h hrj- 2 X lO^m, 29 
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Figure 6: Plots showing the proton flux {(jf + 0*^), the neutron flux (0" + 0*") 
and the total CR flux {(ptotai) versus the energy per particle for a thermal photon 
spectrum with /c^T = 1.8 eV and p^L = 2 x 10^^ eV/cm^. Also shown is the total 
CR flux that is obtained ignoring the escape of neutrons. 

and hence the photon source must be compact (like a pulsar) and in general smaller 
than the typical size of a white dwarf. It has to be mentioned however that in 
non-steady-state situations like in supernova explosions, luminosities larger than 
LiEdd ^-re possible. Another consideration is that if the source has a black-body 
spectrum with temperature T and it is radiating at the Eddington limit, the relation 
Anr^aT^ = Lem implies that T ~ 1.8 x 10^ eV(M/M0)^/4(lO km/ro)^/^ and hence 
for the photon energies to be in the optical/UV region, radii somewhat larger than 
those of neutron stars would be favored (or the luminosity should be smaller than 

LEdd)- 

As a summary, we have considered in more detail the photodisintegration of CR 
nuclei in the scenario in which their interaction with optical and soft UV photons 
around the source is responsible for the steepening of the spectrum and for the 
change in composition above the knee. The more accurate treatment of the nuclear 
processes, in particular the inclusion of the multinucleon emission, implies that 
the required photon column densities are lower (by an order of magnitude) than 
previously obtained. The main difference between this scenario and other (rigidity 
dependent) explanations proposed for the knee is the resulting CR composition, 
which here becomes lighter above ~ 10^^ eV. 
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